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Chapter 1 
Introduction 


High frequency electromagnetic analysis of problems involving complex geometry de- 
scriptions of the scattering bodies has been of great interest in recent years. Spline 
and polynomial type surfaces are commonly used in high frequency modeling of com- 
plex structures such as aircraft, ships, reflectors, etc., and is therefore of interest to 
develop efficient and accurate solutions for the purpose of RCS and antenna pattern 
predictions. Numerical techniques for treating these problems such as Physical Op- 
tics (PO) [1] and the Method of Moments (MM) [2], while simple in concept, still 
require vast amount of computer resources for even intermediate sized targets. The 
use of closed form ray optical solutions, whenever possible, is one way to increase 
efficiency and also provide important physical insight into the various scattering and 
diffraction processes. 

In many instances, high-frequency scattering and diffraction mechanisms can be 
described in terms of ray optical fields that behave locally like plane waves propagat- 
ing along the trajectories of the Geometrical Optics (GO) [3] and the Geometrical 
Theory of Diffraction (GTD) [4]. Although these rather simple ray optical techniques 
can be used to solve many complex practical electromagnetic scattering/radiation 
problems, they fail to describe the more complicated field behavior inside transition 
regions where wave focusing due to higher order phase catastrophes and or changes 
from illuminated to shadow zones occur. Removing the errors within these regions 



would greatly enhance the applicability of GTD as indeed has been shown for its uni- 
form version, the so called UTD [5]. The UTD approach for dealing with transitional 
field behavior involves the use of appropriate transition functions that eliminate the 
errors of GO/GTD solutions, and also provide a smooth connection into the regions 
where these classical ray optical formulations remain valid. 

A common method of analysis for high-frequency scattering and diffraction prob- 
lems involves the use of radiation integrals as well as plane wave integral represen- 
tations for the fields, with the asymptotic approximations for the various scatter- 
ing mechanisms found from the critical point contributions of the integrand. This 
method of analysis when applied to a cubic polynomial surface containing an edge re- 
sults in a stationary phase integral characterized by two stationary phase points that 
are arbitrarily close to one another or to an integration endpoint. Uniform asymp- 
totic evaluation of integrals with such analytical properties involves the incomplete 
Airy function [6] which serves as a canonical integral for the description of transition 
region phenomena associated with composite shadow boundaries. Transition regions 
of such type result from the merging of the reflection shadow boundary associated 
with an edge in the scattering surface, and the smooth caustic of reflected rays aris- 
ing from the confluence of the two stationary phase points near a zero-curvature 
(inflection) point. When the reflection shadow boundary is not in the immediate 
vicinity of the smooth caustic the conventional UTD diffraction coefficient [5] which 
involves the Fresnel integral as a canonical function can be used to effectively de- 
scribe the edge diffracted fields. Furthermore, the ordinary Airy integrals and their 
derivatives are the appropriate canonical functions for the description of the high- 
frequency fields in the neighborhood of the smooth caustic [7]. However, when there 
is a confluence of both reflected and caustic type shadow boundaries, neither the 
Fresnel integral nor the ordinary Airy integrals adequately describe the transition 
region phenomena, and they must be appropriately replaced by the incomplete Airy 
functions. Although these functions are not so easily generated [8] as, for example, 
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the Fresnel integral or the complete Airy functions, a method for their efficient and 
accurate computation was recently developed [9] and allows for the formulation of 
uniform asymptotic solutions that are useful for engineering purposes. 

In this report, an extended UTD solution for the scattering and diffraction from 
perfectly conducting cubic polynomial strips is derived and involves the incomplete 
Airy integrals as canonical functions. The scattering mechanisms involved are de- 
scribed by a new set of uniform reflection, first-order edge and zero-curvature diffrac- 
tion coefficients that remain valid inside the transition regions and also provide 
smooth connection into the regions where the classic ray optical formulations still 
apply. Although additional higher order mechanisms such as creeping waves, double 
edge diffraction, and whispering gallery modes may also exist for the geometries con- 
sidered, they will not be examined in this report. Also, this new solution is universal 
in nature and can be used to effectively describe the scattered fields from flat, strictly 
concave or convex, and concave-convex boundaries containing edges. Numerical re- 
sults obtained using the extended UTD solution showed excellent agreement with 
the method of moments for both polarizations, except for some limited regions in the 
non-specular direction for the TE polarization where the higher order mechanisms 
become significant. 

The outline of this report is as follows: In Chapter 2, a uniform asymptotic 
analysis for the plane wave scattering from a perfectly conducting cubic polynomial 
boundary containing an edge is presented, and in Chapter 3 the extended UTD 
solution for the scattering and diffraction from cubic polynomial strips is formulated. 
In Chapter 4, some indicative numerical results are presented and discussed with 
their accuracy confirmed via comparison with reference moment method results. In 
addition, results obtained using classic UTD are also shown to illustrate the need 
for the new solution. Finally, some concluding remarks and the accomplishments of 
this work are summarized in Chapter 5. 
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It is assumed that all Helds are lime harmonic with time dependence e ?u; ', which 
will be suppressed throughout tliis report. Also, the medium surrounding the scat- 
tered is assumed to be free space. 
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Chapter 2 

Uniform Asymptotic Analysis 

In this chapter, a uniform asymptotic analysis for predicting the scattered fields 
from perfectly conducting cubic polynomial strips is presented. The relevant canon- 
ical geometry for the analysis that follows is a semi-infinite two-dimensional cubic 
polynomial boundary illuminated by a plane wave as shown in Figure 2.1. The scat- 
tering mechanisms to be examined, namely specular reflection, zero-curvature and 
first-order edge diffraction are illustrated in Figure 2.2. As a plane parametric curve 
the boundary C is given by: 

r'(Q') = xx + yy(x) ] a < x < oo (2-1) 

where 

y(x) = a u + a,x + a 2 x 2 + a^x 3 . (2-2) 

The zero-curvature point, x p , is a root of the second derivative of the surface, i.e., 
y"(x v ) = 0 and is given by: 

x (2.3) 

P 3a, ' 

The scattered fields at any point P away from the surface boundary can be 
expressed in terms of the usual radiation integrals over the electric current J g induced 
on the boundary by the incident plane wave as follows: 

E*(P) ~ — [ [s x 3 x J,(Q')]ff, ( , 2) (*r)<tt' (TM case) (2.4) 

v 4 Jc 

ii-lP) ~ -J I /,(<?') x i)H^\kr)di' (TE case) (2.5) 

4 Jc 

5 




Figure 2.1: Canonical geometry for the uniform asymptotic analysis 



Figure 2.2: Scattering mechanisms associated with a cubic polynomial boundary 
containing an edge. 
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where 


r = 
s = 

J,(Q') = 

H^(kr) = 
Z 0 = 
k = 
dl' = 


\Q'P\~P-t' (Q') * 3 (far-zone approximation) , 
x cos 6 + y sin 6 , 

value of J s at any point Q' on the surface, 

cylindrical Hankel function of the second kind of order zero, 

impedance of free space, 

wavenumber of free space, and 

line integration element. 


In the PO approximation to (2.4) and (2.5), the induced current is assumed to be 
given by GO as follows: 


£(<?') 



2n x H'(Q f ) , on the lit portion of the boundary 
0 , on the shadowed portion of the boundary 


( 2 . 6 ) 


in which n is the outward normal to the surface and H'(Q') is the incident magnetic 
field at Q' that under the plane wave incidence assumption is given by: 


| Z^{z x (™ case ) 

[ zeW '^ (TEcase) 


(2.7) 


where 


s l = x cos 9 r + y sin 9 f . 


(2.8) 


Using the assumed current in (2.6) and the large argument form of '(kr) given 
by: 

/ o; 

(2.9) 


H^ 2 \kr) ~ e~’ tr , for kr » 1 


the far-zone scattered field assumes the form: 

Hk e~i kp 

u;(p-,e',6)~T\i r e l ,, h (6',6) 


2tt 


Vp 


( 2 . 10 ) 
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where e< 

,.;,(0',0) is the 

angular field dependence and is given by the stationary phase 

integral 






poo 

«.*(*',*)=/ f,.l ,(T)e**Wix 
J a 

(2.11) 

where 





fs,h{ x ) 


(2.12) 



= f'(Q') • (? + s) = xC{9\ 6) + y(x)S(e\ 6) , 

(2.13) 


c(d\e) 

= cos 6 f + cos 6 , 

(2.14) 


S(0\6) 

= sin 9 f -f sin 9 , 

(2.15) 


dl' 

dx 

= h{x) = \Jl + [^(z)] 2 , and 

(2.16) 


n'(Q') 

iy'( x ) - y 

h(x) 

(2.17) 


The next step in this procedure is the uniform asymptotic evaluation of e,,/,(0',0) 
given in Equation (2.11). 


2.1 Uniform Asymptotic Evaluation of e s ^{0\B) 

The integral in (2.11) can be transformed into a canonical form by first expanding 
the phase function in a Taylor series around the zero-curvature point x p , i.e,, 

<j>{x) = <f>{x v ) + (x - x p )<f>'(x p ) + — ^-0 w (*p) • (2.18) 

Notice that the above expansion is exact since 4>"{x p ) — 0 and <f>( n \x) = 0 for every 
x with n > 3. Next we make the following linear transformation: 

(j>(x) = t(s ) = a + (3s + a 3 /3 , (2.19) 


where 


a = r(0) = <j>(x p ) , and 
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( 2 . 20 ) 

( 2 . 21 ) 



The proper branch for /3 depends on the sign of and <f> (®|j)> thus 




±\ 0 \e-n if rM < 0 f ± foI ^ (Xp) >0 

±|/?| if <j>"'(x t> ) > 0 

(2.22) 


where 

1 


— 


w - M»,)l ^ • 

(2.23) 


Also, alternative expressions for a and (3 may be found directly from Equation (2.19), 

l P 



a = r(si, 2 ) + jj(-/3) 5 

(2.24) 



2 

P = -{^K0)- r ( s '.2)l} 

(2.25) 


where s lt 2 = ±(— /3) l//2 

are the stationary phase points satisfying r'(-S|, 2 ) 

= 0. Thus, 

— 

using Equation (2.19), 

the scattered field angular dependence, e.,,h(0\0), becomes: 

— 

where 

e,,„(<M)= f '°°’ 

''{.a 

(2.26) 



£a — ( a X 1 j) 2 ’ 

(2.27) 



d>X 

G,, h {s) = /,,h( x )^j> and 

(2.28) 

— 

- 

dx 2 1 » 

di “ [r'M. 

(2.29) 


The proper branches for and ^ depend on the sign of (f>"'{x p ), thus 


— 

e. = - 

±l{ " |e ’ 5 if rM <« i ±[0 , <Xrt or 

±ic.i >f rM > 0 

(2.30) 


l£>| = l ffl 

, *'"M i U r, , ,,, \|1 /! 2 S 

2 -I'M"-) < / > ( x p)I ’ 

(2.31) 
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dx 

ds 


- < 


1 

1 d,« 

e 

if 

if 

dr 

ds 



if y(*> 


dx 


2 


ds 




, and (2,32) 


(2.33) 

Next we employ the Chester et al. expansion [10] for the amplitude function in 
Equation (2.26), i.e., 

oo 

G.M = T. KV + 0) m + + 0)"l (2 34) 

m=0 

and since only the leading terms in the asymptotic expansion of (2.26) will be re- 
tained, Equation (2.34) may be written as follows: 

G s ,/i(-s) = a'u' h -f sby' 1 ' -f (s 2 + /?)<7 Si /,(s) 


where 


1 


®u* — 2 1 ) + G^./i^)] ) 




1 

2s, 


[f^s,/i(^i) (^2)] ? and 


g.M = EkiV + r’ + ^V + r 1 ! 


(2.35) 

(2.36) 

(2.37) 

(2.38) 


m = l 


Now, using the far-zone observation approximation and the symmetry of the surface 
near and around s = 0 or x = x p , Equations (2.36) and (2.37) simplify as follows; 

dx I 


a i) h — @sji{ s l, 2 ) — fs,h{% 1 , 2 ) 


ds 


, and 


« = *1,2 


where 


K' h - 0 


dx 

ds 


5=s t|2 


fW ( J 1.2) 


±2 (-/ 3) 1 / 2 

<£"(*1,2) 


(2.39) 

(2.40) 

(2.41) 


and using Equations (2.32), (2.33), (2.39), and (2.41) we have 

K'‘|e- J i if (f>"'{x p ) < 0 


s,h 

(Li i — 


do' 1 1 if^"(*p)>0 


(2.42) 
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where 


i = i/., !.(*i.2)i 


20' 12 


<^"(*1,2) 


= I/*, *(*l, 2) | 




(2.43) 


Then, using Equations (2.35) and (2.40) e,,h(0',0) may be given by the following 
expression: 


' A) = 


e*,h{0 ,8) = e 


J !'■" is + f~‘(s‘ + 0)g,, h {>y HI ” W/: '' 1 is 


(2.44) 


The integrand of the second integral on the right is regular over the entire path of 
integration, thus integration by parts yields 


e.,,i.(0',0) = e* T < 0 > |o£* jT °° e j W' W/:i) ds - ^ 1 ,, h (6)e iW+ ^ /3) 
and using Equations (2.19) and (2.35) it simplifies further 


+ 0(fc- 2 ) 
(2.45) 


..,,(«',#) ~ e'‘" !ni a; i ' 1 [ °° e Ml, ’*’’ l2) ds + e i * T(w 




jfc(/? + ^a) jkr’ita) 


(2 .46) 


where the terms of order higher than 1 / k have been omitted. Next, we let s = /r 1 ^ 

0 = fc- 2 / ;, /3 and in (2.46), thus 

G'At*) 

jkr'{( a ) 

(2.47) 

Finally, the field angular dependence is given by the following expressions depending 
on the sign of <f>"\x v ') and fi(x p ): 


e„ ih (6\0) ~ e jkT( " ) k-' / \ s [) ' h / °° e j( ^ +,,/:! ) dt + e 3 

* ta 


jkT(Za) 


k-'f-'a'* 

IM + C) 


1. 4>"'(x v ) < 0 and </>'(x,,)<0 


In this case, 0 = ae 3 3 , = C a e J 3, a o ,,, = l a o’ , 'l e * 3 an< ^ Equation (2.47) be- 

comes: 


?,, h (e\e) ~ e jH ^h- L *\a?; h \Ai’( 



j f*A a ) 
j( * ~ C 2 ) k 


(2.48) 
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2 . </ ,, ”( r r) 0 and < f > '( T v )> 0 


In this case, (3 = cr, af/‘ = |a, , ) ' , '| and Equation (2.47) becomes: 


e, %h (6\6) ~ e jk ^k-*\aZ h \Ai(*,( a ) + e jU(o) 


fc *K*I , i. /-.*(“) 
* * 


. (2.49) 


[i(<r + C) * *(«) 

Both a and £ a are real quantities with cr <0 if ^ ; (®i>)<® an( ^ Ca<0 if a < x p • Ai(<r, 7 ) is 


the incomplete Airy integral defined by: 

Al(cr, 7 )= f°° e j{ "* +23/3) dz. 

2.2 Total Field Solution 


(2.50) 


Using Equations (2.10), (2.48) and (2.49), the total scattered field is given by: 


1 . <f>'"{x v ) < 0 


u;(p) 


T 

1 e J k, H a 




, e J T fei 

V^ 7 T 


a o h |Ai’(<r,f 0 ) 


k^> |a"/*|ed 


e 3 4 


/*,/.(«) 


j\Z2n{a + Q) y/2n k 


-jkp 


Vp 


with cr (' 0 if tj>\x j,)< 0 . 


(2.51) 


2 . ^'"(x„) > 0 


WM 


jk4>(xp 


t 

_L e W") 


^K'^Aitcr.C.) 

fc5|ay , '|e J T /,./,(a) 






(2.52) 


J j'v / 27t(o’ + Ca) \/ 2 irfc <£'( a ) 

with cr^O if </} r (x f> y^0. 

The in front of the expressions in (2.51) and (2.52) corresponds to the soft or 
hard polarization cases. Notice that the last two terms in Equations (2.51) and (2.52) 
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are purely edge diffraction contributions where the first term in (2.51) and (2.52) 
appears to be a contribution from the zero-curvature point. In fact, three different 
scattering contributions, namely reflection, edge and zero-curvature diffraction are 
implicitly contained in the incomplete Airy integral. It is advantageous to separate 
the scattering mechanisms, extract the appropriate reflection and diffraction coeffi- 
cients, and cast the total solution in a UTD format. In this form, the solution would 
be applicable to more general problems and also provide important physical insight 
into the individual scattering processes. The details of the separation procedure are 
presented in the next chapter. 
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Chapter 3 

Extended UTD Solution 
Formulation 


The total scattered field from the semi-infinite perfectly conducting cubic polynomial 
boundary of Figure 2.2 can be expressed as follows: 

v;( P ) = uhp) + v'M + u%) (3.i) 


where U r z , U c z and U z are the reflected, zero-curvature and edge diffracted field 
components. Although additional higher-order mechanisms such as surface waves, 
double diffraction and whispering gallery modes may also exist, they will not be 
considered in this analysis since the PO approximation to the induced surface cur- 
rents does not account for such effects. Before proceeding with the expressions for 
each field component, the various field contributions that are implicitly contained 
in the incomplete Airy integral are extracted by deforming the original contour of 
integration into appropriate steepest descent paths through the critical points of the 
integrand. 

First, let’s consider the complex plane topology of the incomplete Airy integral 
when the saddle points are real ( cr < 0 ), as depicted in Figure 3.1. There exist three 
cases depending on the location of the endpoint 7 relative to the saddle points z } 
and z^, and Ai(<r, 7 ) may be written as follows: 


Ai (^>7) = / 

J ■> 


00 exp(jV’2) 




/3) dz+ [ ( iz 

Jl™ 
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Figure 3.1: Complex plane topology and contour deformation for the incomplete 
Airy integral when a < 0. 
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(3.2) 


+ / is if 7 < 

^31 


*'3 1 

= r’**'' '*..+* mi'+f e *~+sm dz 

J'y J Ijx | 

if — ( — cr )'/ 2 <7 < (— cr) 1 / 2 , and 

= I'"'" 1 ’'"' c< l '' , ‘*‘' / ' ] dz if 7 > ( - <r) 1/2 . 


(3.3) 

(3.4) 


Using the definitions for the ordinary and incomplete Airy functions found in Ap- 
pendices A and B the integrals in (3.2)— (3,4) may be expressed in terms of these 
functions as follows: 

f»oocxp(i02) 


^ ' ' e ] (<rz+z* /3) fa _ g 2 (o- j7 ) t 

f = 7 r[Ai(<r) — jBi(cr)] , 

•Z/,23 

/oocxpfjV’j) , 

e^^dz = g,(<r, 7 ), 


l 


.11 

oc cxp(ji/’i ) 


/ e i(^+^/ 1 ) = 7 r[ Ai(cr) + jBi(<r)] , and 

^.ii 


e ]("t+‘*/3) dz - g| ((T,7). 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 


Thus, using Equations (3.2)--(3.9), the total scattered field solution in (2.51) and 
(2.52) when <r < 0 (lit side of the caustic) is given by: 


1. < 0 


~ -F (^"^*i|oS- fc | - j’Br (crJlvK-o ')'/ 2 - C-l 

+ c i|(-‘ r )* /? [Ai(<r) + ;'Bi(a-)]u[— (— <r )'^ 2 — £,]} 


+ e 


JW(«0 / ^ I 


e J 4 




o- + Ci 


i) 


-jkn 


Vp 


\/2nk 

2. 0 w (s„) > 0 

tf;(p) ~ T ( V /|e i USW' , ’|{e>**l''le>l<-') , ' , [Ai(»)4jBiW]«[(-<7)'' 2 -(„) 


(3.10) 
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Figure 3.2: Complex plane topology and contour deformation for the incomplete 
Airy integral when a > 0. 

[Ai*(cr) — jB\*(a)}u[—{—<ry ^ 2 — C«] } 


+ e 


Jkd>(a) 


e J 5kz 


| 


\/ 2 n 

\ 

e~ } 5 f*.h{a) 




sfl/k <f>'(a) 


[<? + Q 

e~ jkp 




sfP 


(3.11) 


where 


g 2 (<r,Cc) if C» < -(-^)' /2 

hM) = <1 g 3 (<r,C a ) if -(~<r)'/ 2 < („ < (~*) ,/2 ( 3 - 12 ) 

gl(o’.Cn) ^ Ca > (-< r ) ,/2 

and u(i) is the Heaviside unit step function. The first two terms in Equations (3.10) 
and (3.11) represent uniform reflected field contributions where the last two terms 
represent edge diffraction contributions. 

Next, let’s consider the complex plane topology of ^ST(<r,7) when the saddle points 
are complex ( a > 0), as depicted in Figure 3.2. There exist two cases depending on 
the endpoint location relative to the zero-curvature point (z = 0), and the mcom- 
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plete Airy integral may he expressed in terms of the complete and incomplete Airy 
functions in a fashion similar to the previous case, i.e., 


^T(o , ,7) = g 2 (<r,7) + 27rAi(<r) if 7 < 0, and 
= gi(o‘»7) if 7 > 0 . 


(3.13) 

(3.14) 


Then, using Equations (3.13) and (3.14), the total scattered field of (2.51) and (2.52) 
when a > 0 (dark side of the caustic) is given by: 


1 . < 0 


U;(p) ~ ^(e jk *^V2^e j U L «\a° 0 h |Ai*(ff)ti(-< 0 ) 


+ e 




e j n* 




I ° 


V2nk <P'{a) 


)) 


.*■ + c? 

e -jkp 


+ c i(«rC.+CS/--.) J; (^ Co ) 


Vp 


2 . f'M > 0 


U*(p ) ~ =F |ay /T |Ai(cr)u(— ^ a ) 


+ e 


jk<f>{a) 


e 3 4 k fi 


«u I 


e 

72 


± fM] \\ 

\i vk _ J / 


v + Q 

jkp 


+ c -iK-+Cj/3 J/^^) 




where 


/</(*, 0 = < 


(3.15) 


(3.16) 


(3.17) 


g 2 (°-,Ca) ifCa<0 

gi Ca) otherwise 

In this case the first term in Equations (3.15) and Equation (3.16) represent 
zero- curvature diffraction contribution where the last two terms in (3.15) and (3.16) 
represent edge diffraction contributions. Notice that the zero-curvature diffraction 
contribution is a uniform version of the complex ray (evanescent) field interpretation 
of Ikuno and Felsen [11, 12], and it appropriately reduces to their expressions when 
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<r > 1, The main advantage of our formulation is of course uniformity, and also 
the need for a complex extension of the reflecting boundary and the evaluation of 
surface parameters in complex space is avoided. 

Another observation concerning the edge diffracted field component of the to- 
tal field is that the last term in Equations (3.10), (3.11), (3.15) and (3.16) can 
be recognized as the PO half-plane diffraction coefficient, and using the following 
expressions: 




n°) 


h(a)h'(Q' a ) 


i' 

3 



sin tp' 
sin <p 


h(a)t'(Q' a ) ■ (i* + J) - h(a) C(<p',<p) 


(3.18) 

(3.19) 


it may be written as follows: 


nPO, / , , e->T (sin y>' ; sin <p) 

”' h ^ ^ \Z2nk (cos <p' + cos p) 


(3.20) 


where p' and p are the angles of incidence and observation, respectively, measured 
from the edge-fixed coordinate system. It is well known that the PO diffraction 
coefficient does not satisfy reciprocity nor does it satisfy the local boundary condition 
on the surface for the TM polarization case. These shortcomings produce errors in 
the total scattered field for observation points away from the optical boundaries. 
A rigorous method of providing the necessary corrections to the PO approximation 
involves the introduction of an additional non-uniform induced current near the edge. 
A rather heuristic but simpler approach for improving the PO half-plane diffraction 
coefficient is to introduce a pair of edge correction multiplication terms derived by 
James [13], i.e., 


C s {p',p ) = sin(v?/2) sec(y>72) , and (3.21) 

C h (p',p) = cos(//2) csc(ip/2) (3.22) 


which correct the PO diffraction coefficient so that it yields the exact diffraction 
coefficient outside the optical boundaries. These correction factors reduce to zero 
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Figure 3 . 3 : Geometry for the reflected field from a cubic polynomial strip. 

near the optical boundaries, and therefore they do not afFect the cancellation of 
singularities or the compensation of discontinuities in the reflected field components. 

A final comment about the total field solution in Equations (3.10), (3.11), (3.15), 
and (3.16) is that it remains valid for any value of the cubic polynomial coefficients 
a u , a 1? a 2 , and aj and one can show by carefully taking the appropriate limits that 
it reduces to the half-plane solution when a 2 and a a go to zero, and also reduces to 
the solution for a semi-infinite parabolic screen when a , 3 goes to zero. 

We are now ready to derive expressions for the individual scattered field con- 
tributions from a finite cubic polynomial surface using the expression for the total 
field from the semi-infinite surface given in (3.10), (3.11), (3.15) and (3.16) and by 
expressing the quantities a and £ 0 in terms of the local surface parameters 

relevant to each scattering mechanism. 

3.1 Reflected Field Solution 

The geometry for the reflected field appears in Figure 3.3. The reflected field con- 
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tribution consists of two specular components that only exist in the lit side of the 
caustic. It is given by the following expression: 


U r z {p) ~ u:'(p)u{b -x ri ) + U?{p)u{x Ti - a) 


(3.23) 


where Ul'' 2 (p) are the two specular contributions given by 


/l f 

ui'Hp) ~ 2 ^( 1 ?:,, ,)*'(<?:,.,) - 

R,,{Qr ) is the principal radius of curvature at the reflection point given by 

R (O’ ) = 

and * r , 2 are the specular points found by 


y/P. 


(3.24) 


(3.25) 


X r>,l ~ 


— 02 dh Y a 2 — 
2d} 


(3.26) 


where 


a, = i • (s' + s r ) + a,# • (s' + s r ) = (7(0', 0) + aiS(0', 0) , (3.27) 

a 2 = 2a 2 y-{s' + s r ) = 2a 2 S(8',e), and (3.28) 

a,i = 3«3»-(i i + 0 = 3a 3 5(^,«). (3-29) 

The quantity 7 Zs,h(Q' ri ) 1S the uniform acoustic soft or hard reflection coefficient that 
remains valid as Q' r . -» Q' p and is given by 


n.M' ri ) = R,,h 


TrM if R a {Q' ri )> o 
t;K) ifi? g ((?' rj )<o 


(3.30) 


where 


R Sy h 

T r (* r ) 

<T r 


== fpl , for a perfectly conducting boundary, 

= [Ai(-o-, ) - J*Bi( — cr r )] , and 


h[f\Q' ri )-f'(Q' p )).(s’ + n 


(3.31) 

(3.32) 

(3.33) 
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dark side jj / lit aide 

caustic 

Figure 3.4: Geometry for the zero-curvature diffracted field from a cubic polynomial 
strip. 


The quantity T r (a r ) is the appropriate transition function that cancels the GO field 
singularity at the caustic. Also, when the specular point is far removed from the zero- 
curvature point (observation in the deep lit side of the caustic), <r r >> 1 and using 
the large negative argument forms of the ordinary Airy functions (see Appendix A) 
we have that T, (<r r 1) ~ 1 and — ► R^h so that the extended UTD reflected 

field solution of Equation (3.24) reduces to the usual GO expression for the reflected 
field. 


3.2 Zero-Curvature Diffracted Field Solution 

The geometry for the zero curvature diffracted field is shown in Figure 3.4. The 
zero-curvature diffracted field contribution exists only in the dark side of the caustic 
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ami compensates the reflected field discontinuity. Tt is given by 

u;(p) ~ ^ «(*„ - «)„((, - .„) (3.34) 

VP 

where V' sh {Q' v ) is the scalar zero-curvature diffraction coefficient and is given by 






i 


I r .(v ?, ^)Ai(<r c ) 


where 


W,v) 


, , _ j/gV^ + C' 2 (y , ,yj 


I c jC(^',^) - c 2 S(y> , ,y>)| | / :j 


cos 


1 ^) 


(3.35) 


(3.36) 


with ci = 


and c 2 


1 


h{x p ) — ^-A(«,) 

is the diffraction coefficient angular dependance, and 


Jfes 


y"'(* P ) 


|CV»¥>)I 


(3.37) 


c 2 |ci C(^', v?) - c 2 S{if <p )\ '/ 3 ' 

Notice that the zero-curvature diffraction contribution is significant near the caustic 
and decays exponentially away from the caustic as Ai(tr c > 1) — > 0. 

3.3 First-Order Edge Diffracted Field Solution 

The geometry for the edge diffracted field is illustrated in Figure 3.5. The edge 
diffracted field contribution consists of two components (one from each of the two 
edges) and is given by 

U?(p) = U^(p) + Uf\p) (3.38) 

where 


U<»(p) ~ e ifc, " (<? -.‘ ) i ''ir )/) (Q; ii ,)e ifrr '( t? M) i 


d e 


-jkp 


sTp 


(3.39) 


an( I 'D'*,h(Qa,b) i s the scalar uniform edge diffraction coefficient that remains valid as 
Q'„,b ~ > Q',,- It is comprised of two parts, i.e., 


KM'.t) = d'IW:*) + KhiQ'ot) 


(3.40) 
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R.S.B. 


Figure 3.5: Geometry for the edge diffracted field from a cubic polynomial strip, 
where D hp h (Q ' a< ,,) is the half- plane diffraction coefficient given by 


_ e -if f f'f'nl,- fa,b\ ( ‘Pn.b + ¥>«,b\ 

D - 5«W = H 1^-i— J * s ' c ) 

and D r , h {Q' a b) is a term ^ at provides curvature correction and is given by 


(3.41) 


, , % (talkie# 

D's,h(Qa,b) — "F ^2^ 


y w (a,6) 


■^d(^a,6’ < P a &) 


U<r*J,,U) if y’"(a,b)y ■ (s’ + s J ) < 0 

T;^mU) if y"'{ a >b)y ■ (s' + a 1 ') > 0 

(O /10\ 


where 


ThMC\) 

T d (v,-\(\) 

^a,b 


+ ^Wajbi W*>) 

\dlC(<p' a<b ,<Pa,b) ~ ^25'(v?a.6)y ? a.fc)l </3 


COS 


( 4 


.6 V’a.b 


+ e 


cr+ICI 2 
-T d >,|C|), 


= fc5 


sgn 


;h<,W 6 ^’ ICI) if ' + IC, ^° 

gt(^ICI) if o- + Id 2 > 0 

5 |(f;)G(yga y>o,b) + d 4 S{fp' a ji t ,(Pa,b)\_ _ 

|d| b , y’n.b) — d-2S{tp' a ji,,<Pa, 

X • (a 1 + i rf ) + y'(^/))y • (^’ + f^) 
iTMHl- + a rf ) . 


y"'{a,b) 


(3.43) 

(3.44) 

(3.45) 


(3.46) 
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da,b 

d t 


di 


ds 

d.x 


\dlC(<p’ ab ,<p a ,b) - diSfaajjftpafi)^ 5 a i b< x p 1 


±1 ) 

■ y'{a,t>) 

h{a,b) ’ 

1 

h(a, b) ’ 

1 + y'ja^y^x,,) 

h(a, b) 

y'{a,b) ~ j/faQ 

h(a, 6) 


and 


(3.47) 

(3.48) 

(3.49) 

(3.50) 

(3.51) 

(3.52) 


Notice that the curvature correction part of the uniform edge diffraction coefficient in 
Equation (3.40) provides significant contribution near the reflection shadow bound- 
ary (RSB) where a + ( 2 % 0. In fact both the half-plane and curvature components 
in (3.40) become singular near the RSB, however these singularities cancel each 
other out and the total diffraction coefficient remains finite. Away from the RSB 
( 0 - -(- £ 2 0) the curvature correction component reduces to zero as the diffraction 

transition function Tj(cr,() — > 0, and thus the edge diffraction solution reduces to 
the half-plane solution as it should. 
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Chapter 4 

Numerical Results and Discussion 


In this chapter, some numerical results for the scattered fields from a general cubic 
polynomial strip are presented. The scattering geometry and the relevant parameters 
are illustrated in Figure 4.1. The accuracy of the extended UTD solution is verified 
via comparison with method of moments results. Also, some results obtained using 
classic UTD are shown and illustrate the need for the new solution. 

The first example considered is a cubic polynomial strip with a 0 = 2.0A, a { = 0.5, 
a 2 = 0. 1 A “ 1 ? a* — 0.1A -2 , a = — 1.5A, and b — 1.5A. In this case both edges 
are far removed from the zero-curvature point and thus the two reflection shadow 
boundaries and the caustic of the reflected rays are clearly distinct. Figure 4.2 shows 
a plot of the magnitude of the various field contributions to the total field for the 
TM polarization case and an angle of incidence O' = —45°, Notice that the reflected 
field component exhibits a total of three discontinuities. The first discontinuity 
occurs across the caustic of the reflected rays at 6 m —82° and is compensated 
by the zero-curvature diffracted field. The second discontinuity occurs at the RSB 
associated with the edge Q a at 0 % —50° and is compensated by the edge diffracted 
field from Q a . Similarly the third discontinuity occurs at the RSB associated with 
the edge Q\> at 6 « —23° and is compensated by the edge diffracted field from Qi,. 
Also notice that the two edge diffraction terms become singular near the incidence 
shadow boundary (ISB) at 9 — 135°; however, they combine to give a finite result. 
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igure 4.1: Scattering geometry and relevant parameters for the numerical results 




Figure 4.3 shows a plot for the total scattered field in terms of histatic erlio width 
for the TM polarization. The extended UTD result shows excellent agreement with 
the method of moments where classic UTD gives an erroneous result near the caustic 
at ^ si —82°, and in both the lit and dark sides. This is expected since the classic 
UTD formulation does not contain zero-curvature diffraction information in the dark 
side and also uses the non-uniform GO expression for the reflected field contribution 
in the lit side. Figure 4.4 shows results for the bistatic echo width for the TE 
polarization. Again the extended UTt) result shows good agreement with method of 
moments. The discrepancies for observation directions near grazing are attributed 
to higher order mechanisms missing from the total field. These higher order effects 
such as edge exited surface rays and whispering gallery modes are stronger for this 
polarization since the grazing fields do not vanish on the boundary as is the case for 
the TM polarization. The failure of the classic UTD near the caustic is again clearly 
illustrated. Figures 4.5 and 4.6 show plots of the monostatic echo width for the 
TM and TE polarizations, respectively. The extended UTD solution gives accurate 
results for the monostatic case also, except for the regions where the higher order 
mechanisms become significant. Contrary to the classic UTD result, the extended 
UTD solution remains finite and continuous across the caustics. 

For the second example we consider a cubic polynomial strip with a {) = 2.0A, 
a \ = 0-5, a 2 = 0.1A~ 1 , <13 = 0.1A~ 2 , a = — 0.33A, and b = 1.5A. In this case, the edge 
Qn coincides with the zero-curvature point Q r and thus the RSB and the caustic of 
the reflected rays coalesce to form a composite shadow boundary. Figure 4.7 shows 
a plot of the bistatic echo width for the TM polarization and an angle of incidence 
8' = —45°. The classic UTD result exhibits a singularity at the RSB associated with 
edge Q n , where the extended UTD result remains finite and is in excellent agreement 
with the reference solution. Figure 4.8 shows a plot of the monostatic echo width for 
the same geometry. Again the non-uniformity of the classic UTD solution is clearly 


28 



evidenl where tlie extended UTD solution remains valid across flip RSD and reduces 
to the classic UTD result away from the transition region. 

The remaining two examples illustrate the universal nature of the extended UTD 
solution. For the third example we consider a cubic polynomial strip with a tl = 2 . 0 A, 
a, = 0.5, a 2 = 0. 1 A~ 1 , a 3 = 0 . 0 A“ 2 , a = — 0.33A, and b = 1 . 5 A. This of course 
corresponds to a parabolic screen and is well known that classic UTD gives accurate 
results. In this case the zero-curvature point theoretically moves to negative infinity. 
Figure 4.9 shows a plot of the monostatic echo width for the TM polarization and 
the extended UTD solution remains valid and shows excellent agreement with the 
reference solution. 

For the final example we consider a flat strip by letting both the quadratic and 
cubic coefficients go to zero. Figure 4.10 shows a plot of the bistatic echo width for 
the TM polarization and an angle of incidence 9 ( — —45°. The extended UTD solu- 
tion remains valid for this special case also, and clearly demonstrates its flexibility 
for treating general surfaces that are highly curved, slightly curved, or completely 
flat. 
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Field Magnitude 



Angle of Observation Q (degrees) 


Figure 4.2: Scattered field contributions (TM polarization case) from a cubic polyno- 
mial strip with a u = 2.0A, ci| = 0.5, a? = 0.1A -1 , aj = 0.1A 2 , a — — 1.5A, b = 1.5A 
and angle of incidence 6' = —45°. 




Angle of Observation 6 (degrees) 

Figure 4.3: Bistatic echo width (TM polarization case) from a cubic polynomial strip 
with o„ = 2.0A, a, = 0.5, a 2 - O.lA -1 , a 3 = 0.1A" 2 , a = -1.5A, b = 1.5A, and angle 
of incidence 6' = —45°. 
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Figure 4.5: Monostatic echo width (TM polarization case) from a cubic polynomial 
strip with a 0 = 2.0A, a) = 0.5, a 2 = 0.1 A - 1 , a 3 = 0.1A -2 , a = — 1.5A, and 6 = 1.5A. 
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Figure 4.7: Bistatic echo width (TM polarization case) from a cubic polynomial strip 
with a , i = 2.0A, ai = 0.5, a 2 = O.lA -1 , a.j = 0.1A -2 , a = — 0.33A, b = 1.5A, and angle 
of incidence 9' — —45°. 
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Figure 4.9: Monostatic echo width (TM polar: 
strip with a 0 = 2.0A, a\ = 0.5, a<i — 6.1A“* 
— 0.33A, and b — 1.5A. 





Angle of Observation 8 (degrees) 


Figure 4.10: Bistatic echo width (TM polarization case) from a cubic polynomial 
strip with a t , = 2.0A, a,\ — 0.5, a 2 = 0.0A -1 , a 3 = 0.1A 2 (flat strip), a = — 0.33A, 
b = 1.5A, and angle of incidence 6' = —45°. 




Chapter 5 

Summary and Conclusions 


We have presented an extended UTD solution for the scattering and diffraction 
from cubic polynomial boundaries containing edges. This new solution involves the 
incomplete Airy integrals as canonical functions and effectively describes the tran- 
sitional field behavior associated with composite shadow boundaries and caustics of 
the reflected rays. The total solution is presented in a ray optical format by deriving 
the appropriate uniform reflection, zero-curvature and edge diffraction coefficients 
that remain valid inside the transition regions, and also provide smooth connection 
into the regions where the classic ray optical formulations remain valid. 

It was shown by comparison with a reference method of moments solution that 
the extended UTD TM polarization solution yields excellent results. The TE polar- 
ization solution also showed good agreement with the reference solution although it 
would benefit from the inclusion of higher order effects such as edge excited creep- 
ing waves, double edge diffraction, and whispering gallery modes in certain regions 
where grazing fields exist. The failure of classic UTD solution to describe the scat- 
tered fields near caustics of the reflected rays and composite shadow boundaries was 
also clearly illustrated. The universal nature of the extended UTD solution was 
demonstrated by considering examples of a parabolic screen and a flat strip. In both 
special cases, the new solution gave excellent results. Therefore, the extended UTD 
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solution ran lie used to effect i vely describe the scattered fields from general curved 
surfaces ranging from strictly concave or convex, concave-convex, or completely flat. 
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Appendix A 

Ordinary Airy Functions 


The ordinary Airy functions satisfy the following differential equation: 

y"(tr) - o-y(o-) = 0 (A.l) 

which has two independent solutions, Ai(er) and Bi(<r). In integral form they are 
given by 

Ai(cr) = [ e j ^ z+z3/:i) dz (A.2) 

27 r J jj’2\ 

Bi(er) = -L [ e j ^ z+z3/3) dz (A.3) 

where the contours of integration Z 2 i, L 23 and Lj 1 are shown in Figure A.l. For 
small arguments they can be computed using their ascending series form [14], i.e., 

Ai(<r ) = C|/(<r) - c 2 g(cj) (A.4) 

Bi(cr ) = V3[c l f(a) + c 2 g{cr)] (A. 5) 

where 


f(a) = 

, 1 3 1-4 G 1-4-7 9 

i+ r + 6 '*+ 9 * +'•• 

(A. 6 ) 

9{°) = 

(7+ 2 < 7 ' + 2 ' 5 (r ' + 2 ' 5 ' 8 <r , "+--- 
4 7 10 

(A.7) 

c\ - 

Ai( 0 ) = Bi(0)/\/3 = 3- 2/3 /r(2/3) 

(A.8) 

c 2 = 

Ai'( 0 ) = Bi'(0)/\/3 = 3- |/3 /r(l/3) . 

(A.9) 
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Figure A.l: Contours of integration for the ordinary Airy Functions. 
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Tlieir large argument forms are given by 


Ai(cr) ~ 

•'V'Vl' 1 " (kg 

<r| < 7r) 

(A.10) 

Ai( — <r) ~ 

I'-'I'v-'l' sin (|' 7 ' /2 + 

(|arg(r| < 2 tt/3) 

(A. 11) 

Bi(c) ~ 

7r _l/,2 cr -1 ^ ,, e5' T * /3 (|arg<r| 

< »/3) 

(A. 12) 

Bi( — a) ~ 


(|arg cr| < 2 tt/3) . 

(A.13) 


Complete asymptotic expansions may be found in [14]. 
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Appendix B 

Incomplete Airy Functions 


The incomplete Airy functions satisfy the parabolic differential equation applied by 
Fock [15] and others to the study of fields near the surface of convex diffracting 
bodies, i.e., 

r xi ^ i 

y(cr, 7 ) = 0 (B.l) 


_ _ .d_ 
dcr 2 ^ d~j 


which has three independent solutions, g^o", 7 ), i — 1,2,3. In integral form they are 
given by 

g,((T, 7 ) = /° CCXP(JV,) e^ 2+t3/3) dz ; i = 1,2,3 (B. 2 ) 

A 

where the contours of integration are shown in Figure B.l. The functions g 2 and g, 
can be obtained from g, and the ordinary Airy functions as follows: 


g^o^) = — 27rAi(cr) , and 

g:i(^ 7 ) = gi(^> 7 ) — ?r[Ai(<T) + iBi(cr)] 


(B.3) 

(B.4) 


For small values of < 7 , gj can be computed using its ascending series form [9], i.e., 

oc 

gi( <r >7) = X! a »(7)< T " ( D - 5 ) 

where 


n=0 


00(7) = 6,(0, 7) = e J>/6 3“ 2 ^ :1 r(l/3, — jV*/3), 
a, ( 7 ) = A gl ( 0 , 7 ) = -e-^ ,i 3- ,/3 r(2/3,-jV/3), 


(B- 6 ) 

(B.7) 
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Figure B.l: Contours of integration for the incomplete Airy Functions. 
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p 

cr 

II 

CT 

(n.s) 

, , ^>1-3(7) + j a u-Al) 

<M 7 ) - > 71 - 3 

n[n — 1 ) 

(B. 9 ) 


and T(x, y) is the incomplete Gamma function. The asymptotic forms of g, are given 

by 


6i(^7) 


—L-e^^/V (7 > |<r| ,/2 ) 

cr + 7 2 


(B.10) 






y/ne 3 " u(—T]) — 
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F'W). 


>J'r 


2i7 

7(<n+-> 3 / 3 ) 


where 


2j [0- + 7 2 cr ) 1 /^ 

77 = ± 0-7 + 7 3 /3 + ^( _cr ) 3/2 ] » 7 <(~ <T ) ) 


(O’ < -1) 

(B.ll) 

1/2 

(B . 12) 


and F(x) is the Kouyoumjian-Pathak transition function [5]. Complete asymptotic 
expansions for g { along with some representative plots may be found in [9], 
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